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RESUMEN 

Los estudios del Gas Ionizado Difuso (DIG), hasta el momento se han hecho 
sin un consenso del criterio estricto para diferenciar entre el DIG y las regiones 
H II. En este trabajo recopilamos las mediciones de las h'neas en emision de 29 
galaxias disponibles en la literatura, creando la primera base de datos del DIG (DI- 
GEDA). Haciendo uso de esta base, analizamos las propiedades globales del DIG 
a partir de los cocientes de hneas [N II]A6583/Ha, [O I]A6300/Ha, [O III]A5007/H/3 
y [S II]A6716/Ha, asf como la medida de emision de Ha. Este analisis nos per- 
mitio concluir que el cociente [N II] /Ha es un criterio general para diferenciar si 
una region en emision es DIG o region H II, mientras que la EM(Ha) es una can- 
tidad util unicamente cuando consideramos galaxias por separado. Finalmcntc, 
encontramos que las regiones clasificadas como DIG en galaxias Irr parecen tener 
un comportamiento mas cercano al de las regiones H II que al DIG en galaxias 
espirales. 

ABSTRACT 

Studies of the Diffuse Ionized Gas (DIG) have progressed without providing 
so far any strict criterion to distinguish DIGs from H II regions. In this work, 
we compile the emission line measurements of 29 galaxies that are available in 
the scientific literature, thereby setting up the first DIG database (DIGEDA). 
Making use of this database, we proceed to analyze the global properties of the 
DIG using the [N II]A6583/Ha, [O I]A6300/Ha, [O III]A5007/H/3 and [S II]A6716/Ha 
lines ratios, including the Ha emission measure. This analysis leads us to conclude 
that the [Nil] /Ha ratio provides an objective criterion for distinguishing whether 
an emission region is a DIG or an H II region, while the EM(Ha) is a useful quantity 
only when the galaxies are considered individually. Finally, we find that the emission 
regions of Irr galaxies classified as DIG in the literature appear in fact to be much 
more similar to HII regions than to the DIGs of spiral galaxies. 

Key Words: astronomical data bases: miscellaneous — galaxies: general — ISM: 
HII regions — ISM: lines and bands — line: formation 



1. INTRODUCTION 
From observations of free-free continuum ab- 
sorption, Hoyle & Ellis (1963) inferred the exis- 
tence of a new gas component of the interstellar 
medium (ISM). They proposed that this new com- 
ponent was warm (T e ~ I0 4 K), of low density 



10- 



cm 



~ 3 ), ionized and that it surrounded 



the Milky Way's disk. The integrated luminosity 
of this hypothetical medium was comparable to that 
expected from photoionization by O and B stars from 
the disk. 

The existence of this gas component was later 
confirmed by the signal dispersion from pulsars and 



by the detection of faint optical emission lines from 
the galactic diffuse interstellar medium (Reynolds, 
Scherb, & Roesler 1973). Deep images of exter- 
nal galaxies such as NGC 891 show a similar warm 
gas component (Dettmar 1990; Rand 1998). This 
component has since been found in several other 
galaxies (e.g., Bland-Hawthorn, Sokolowski, & Cecil 
1991a,b; Veilleux, Cecil, & Bland-Hawthorn 1995; 
Greenawalt, Walterbos, & Braun 1997; Martin & 
Kennicutt 1997; Wang, Heckman, & Lehnert 1997). 
Today, it carries various names, such as WIM (Warm 
Ionized Medium), DIM (Diffuse Ionized Medium), 
DIG (Diffuse Ionized Gas) or eDIG (extraplanar 
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DIG), and "Reynold's Layer" (in the case of the 
Galaxy). Hereafter, we will use the generic names 
DIG and eDIG. 

Despite the fact that the DIG in our Galaxy 
accounts for ~ 90% of the ionized medium's mass 
(Reynolds 1991), the nature of this low density 
plasma is far from being well understood. The 
first problem encountered in DIG studies is how to 
formally distinguish DIGs from H II regions, which 
similarly consist of a warm and ionized gas phase. 
Several methods have been proposed in the litera- 
ture. Examples are the Ha equivalent width (Bland- 
Hawthorn et al. 1991a), the Ha emission measureQ, 
EM (e.g., Walterbos & Braun 1994), or the sur- 
face brightness (e.g., Ferguson et al. 1996). Nev- 
ertheless, none of the methods has so far been ac- 
cepted as standard, mainly because the density and 
size of the emitting regions are usually unknown. 
To illustrate the above problem, an example is pro- 
vided by the galaxy M31, for which Walterbos & 
Braun (1994) adopted the EM criterion to distin- 
guish DIGs from H II regions. They found that close 
to low star formation regions, the DIG always had 
an EM < 50pccm~ 6 , while near higher star forma- 
tion regions, the latter had an EM near lOOpccm -6 . 
Therefore, even though the EM is proportional to the 
electron density and should thereby reflect how "dif- 
fuse" an emission region is, this criterion was shown 
to vary with position in this galaxy. 

Another problem in the study of the DIG is the 
difficulty in observing it, mainly because of its low 
surface brightness. As a result, the observations in 
the literature are sparse and only a few line ratios are 
typically reported. Furthermore, systematical com- 
parative studies of the general behavior of the EMs 
or of line ratios among different objects are still miss- 
ing. For example, even though a majority of authors 
agree on the existence of an anti-correlation between 
the [N II] A6583/Ha or [S II] A6716/Ha ratios with the 
EM (these ratios increase when the EM decreases, 
e.g., Bland-Hawthorn et al. 1991a; Rand 1998), very 
few studies were done with more than two galaxies 
(e.g., Wang et al. 1997; Zurita, Rozas, & Beckman 
2000; Tiillmann & Dettmar 2000; Miller & Veilleux 
2003). In this work, we will verify for the first time 
whether such an anticorrelation is a general feature 
of the DIG or not. Another interesting aspect is the 
apparent confusion concerning the [O III]A5007/H/3 
ratio in the literature, for which some authors re- 
port a constant behavior with EM (e.g., Wang et al. 
1997, at least for some galaxies), while others report 

lr The EM in units of peem -6 is defined as EM = 
J n e n H + ds «s J rig ds. 



an increase when the EM decreases (e.g., Rand 1998; 
Tiillmann & Dettmar 2000). 

The spectral characteristics (different from those 
of H II regions) combined with the large spatial ex- 
tension of the DIG (typically \z\ « 2 kpc), are not 
well reproduced by models that consider only pho- 
toionization from O-B stars (e.g., Rand, Wood, & 
Benjamin 2009), even though it could provide the 
entire UV photon flux required to mantain the DIG 
ionized (e.g., Reynolds 1990; Zurita et al. 2000). 
As a result, several authors have suggested that 
the DIG could be excited by extra ionization (or 
heating) sources like: shocks (e.g., Sivan, Stasinska, 
& Lequeux 1986), mixing layers (e.g., Binctte ct 
al. 2009), old hot stars (e.g., Sokolowski & Bland- 
Hawthorn 1991), and decaying dark matter (e.g., 
Dettmar & Schultz 1992). In the study of the Seyfert 
2 galaxy NGC 1068, Bland-Hawthorn et al. (1991b) 
proposed photoionization by the AGN power law. 
None of these models is completely successful in ex- 
plaining the line ratios, spatial extensions and lumi- 
nosities of the DIGs in spiral galaxies. For a com- 
plete review of the posibles DIG ionizing or heat- 
ing sources, see Bland-Hawthorn, Freeman, & Quinn 
(1997). 

With the previous arguments in mind, and to 
correct for the lack of comparative studies in the 
field of diffuse emission, we have compiled a com- 
prehensive database of emission line ratios or fluxes, 
from DIGs or HII regions, published in the litera- 
ture. This database is hereafter called the Diffuse 
Ionized Gas Emission Database (DIGEDA) which 
will soon be available at the VizieR website of CDS 
(Centre de Donnees Astronomiques de Strasbourg, 
http : // vizier . u- strasbg . f r ) . 

The structure of the article is the following: after 
describing DIGEDA in § [21 we proceed with the 
analysis and interpretation of the data in § [3l A 
discussion follows in § [4] and a brief conclusion is 
given in § [5j 

2. THE DIG DATABASE 

An extensive research on DIG observations in the 
literature led us to build a comprehensive database 
comprising DIG and H II region spectroscopic ob- 
servations of 29 different galaxies (25 spiral galaxies 
and 4 irregulars) . This survey contains galaxies with 
a significant spread in star formation rates, Ha lu- 
minosities, distances, disk inclinations, slit positions 
and slit orientations. 

2.1. Internal data structure of DIGEDA 

Table Q] lists the 17 bibliographic references used 
in DIGEDA, which amounts to a total of 1061 line 
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TABLE 1 



BIBLIOGRAPHICAL REFERENCES USED FOR THE COMPILATION OF DIGEDA 



R.ef_ID 


Author 


Galaxy 


Ollt 


Observations 


No. regions 


1 


Bcnvcnuti, D'Odorico, & Pcimbcrt (1976) 


M33 


Intcrarm 


[O in] A4959, [O in], [N ir] 


62 


2 


Golla, Dcttmar, & Domgorgcn (1996) 


NGC 4631 


X 


[Nii], [Sn] 


39 








II 


[Nil], [So] 


80 


3 


Grccnawalt ct al. (1997) 


M31 


DIG, Hn 


[On], [Om], [Nn], [Sn] + , EM 


18 


4 


Martin & Kcnnicutt (1997) 


NGC 1569 


DIG 


[Om], He i, [Oi], [Nn], [Sn] + 


25 


5 


Wang et al. (1997) 


M51 


Intcrarm 


[Om], [Nn], [Sn] + 


3 






M101 


Intcrarm 


[Ont], [Nn], [Sn] + 


6 






NGC 2403 


Interarm 


[Om], [Nn], [Sn] + 


5 






NGC 4395 


Intcrarm 


[Om], [Nn], [Sn] + 


3 


6 


Rand (1998) 


NGC 891 


II 


[Oi], [No], [Sn] 


45 








± 


[Om], [Oi], [No], [S n] 


in 


7 


Galarza, Walterbos, & Braun (1999) 


M31 


DIG, SBRs 


[On], [Om], [Nn], [Sn] + , EM 


30 


8 


rioopcs, walterbos, &c nana ^lyyyj 


Mpp A (? Q 1 

IN Kji^ 40oi 


_L 


[\J 1 1 1 J , IS IlJ-r, HjM 


31 


9 


Tiillmann & Dcttmar (2000) 


NGC 1963 


_L 


[Om], He i, [Oi], [No], [Sii], EM 


40 






NGC 3044 


_L 


[Om], He i, [Oi], [Nii], [Sii], EM 


47 






NGC 4402 


± 


[N n], [S n], EM 


24 






IN GO 4bo4 


X 


[N I Ij , [b II J , JjjM 


40 


10 


Collins & Rand (2001) 


NGC 4302 


_L 


[Nn], [Sn] + 


13 






NGC 5775 


_L 


[Oi], [Om], [Nn], [Sn] + 


31 






ULrU 102OO 




[vJIJ, [UIIIJ, [1NIIJ, [o IIJ + 


19 


12 


Otto, Gallagher, & Reynolds (2002) 


NGC 5775 


_L 


[On], [Om], [Nn], [Sn] + 


38 






NGC 4634 


_L 


[On], [Om], [Nn], [Sn] + 


34 








II 


[On], [Om], [No], [S n] + 


32 






NGC 4631 


II 


[On], [Om], [No], [Sn] + 


34 






NGC 3079 


II 


[On], [Om], [Nn], [Sn] 


53 






NGC 891 


II 


[Oil], [Nn], [Sn] + 


50 


13 


Hoopes & Walterbos (2003) 


M33 


Interarm 


[Om], He i, [No], [So], EM 


37 






M51 


Interarm 


[Om], [No], [Sn], EM 


26 






M81 


Interarm 


[Om], [No], [So] 


2 


14 


Miller & Veilleux (2003) 


UGC 2092 


_L 


[Nil], [So] 


4 






UGC 3326 


_L 


[Nii], [Sn] 


6 










1 1 1 J , nc I, [UIJ, [1NIIJ, [o IIJ 


13 






NGC 2820 


J_ 


[O m], He i, [Oi], [Nii], [S ii] 


13 






NGC 3628 


_L 


[Om], [No], [So] 


13 






NGC 4013 


_L 


[Om], He i, [Oi], [Nii], [Sii] 


9 






NGC 4217 


± 


[Om], [Oi], [Nii], [Sii] 


17 






NGC 4302 


± 


[Om], He i, [Oi], [Nii], [Sii] 


9 






NGC 5777 


± 


[Om], [No], [So] 


6 


IS 


Hidalgo-Gamez (2006) 


ESO245-G05 


DIG, Hn 


[On], [Om], [Nn], [Sn] 


7 






Gr8 


DIG, Hn 


[Om], He i, [Nn], [Sn] 


19 


16 


Voges & Walterbos (2006) 


M33 


DIG 


[Om], [Oi], [Nii], [Sii] 


7 


17 


Voges (2006) 


NGC 891 


II 


[Om], [Oi], [Nii], [Sn], EM 


10 








± 


[Om], [Oi], [Nii], [Sn], EM 


10 


18 


Hidalgo-Gamez (2007) 


DDQ53 


DIG, Hn 


[Om], [Sn] 


26 



a These numbers also serve as bibliographical entries in DIGEDA. 

b Slit orientation in each galaxy. The symbols and descriptors are described in § 12.11 

c The line wavelengths are the following: [Oll]A3727, [OIII]A5007, HelA5876, [Ol]A6300, [NII]A6583, [SII]A6716 and 
[S II]A6716 + A6731. 



measurement data sets. Table 1 also summarizes the bibliographic references, Column 3 lists the names of 

main characteristics of the observations contained in the galaxies studied in each reference, Column 4 is 

each reference. The table is structured as follows: the slit orientation with respect to the galactic plane 

Column 1 attributes a unique identification number in the case of edge-on galaxies: parallel (||) or per- 

to each reference, which is later used as a reference pendicular (_L), while for the case of face-on galax- 

label in DIGEDA (ReLID), Column 2 contains the ies, a short descriptor of the region is used instead: 
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TABLE 2 
COLUMN HEADERS IN DIGEDA a 



Column 





1 


2 


3 


4 


Information 


CDS_iJJ 


Position 


TO TTl \ Q707 
[{J 11J AO 1 Z 1 


TO TTTl \ A QfiQ 
[KJ 111J A4000 


rip 


Column 


5 


6 


7 


8 


9 


Information 


[OIII] A4959 


[O III] A5007 


He I A5876 


[OI] A6300 


[Nil] A6548 


Column 


10 


11 


12 


13 


14 


Information 


Ha 


[Nil] A6583 


[S II] A6716 


[SII] A6731 


[S II] A6716+6731 


Column 


15 


16 


17 


18 


19 


Information 


T c (10 4 K) 


[S II] A6716/6731 


[O III] A5007/4959 


Ha/H/3 


EM 


Column 


20 


21 


22 


23 


24 


Information 


RefJD 


Morphology 


Slit 


RcgionJD 


GalJD 



"See tj b.ll for further details about the information contained in the various entries. 



Interarm, DIG (Diffuse Ionized Gas), HII (HII re- 
gion), SBR (Super Bubble Remnant). In Column 5, 
we detail the set of emission lines reported for each 
galaxy while in Column 6 we give the total number 
of areas or positional observed within each galaxy. 
There are two factors that may affect the accuracy 
of the line ratios: dust reddening and the presence of 
underlying absorption lines. We have not performed 
any reddening correction although we did use the 
reddening-corrected ratios when made available by 
the authors. In any case, all line ratios considered 
in this Paper involve lines that are relatively close 
in wavelengths and are therefore not very sensitive 
to reddening. As for the presence of underlying ab- 
sorption lines, this is a concern only for Ha and H/3. 
The equivalent width of either absorption line is es- 
timated to be ~ — 2 A (Antonio Peimbert, private 
communication). Since the emission line equivalent 
widths are not available from the works cited in Ta- 
ble HJ we could not correct the ratios for this effect 
nor evaluate precisely its impact. Further details 
about the data reduction performed by each author 
can be found in the corresponding bibliographical 
entry listed in Table Q] Hereafter, we will consider 
the symbols [Nil], [OI], [OIII], [SII], and [SII]+, as 
referring to the [NII]A6583, [Ol]A6300, [OIII]A5007, 
[S II]A6716 and [S II]A6716 + A6731 emission lines, re- 
spectively. 

The 1061 observations obtained from these ref- 
erences were extracted by digitalization of published 
figures or tables. The data were subsequently nor- 
malized and incorporated in DIGEDA. This re- 
sulted in a table of 25 comma-separated columns 
(csv format) containing 1061 data lines (or records). 

2 Throughout the paper, we also refer to these positions as 
"regions" . 



Missing entries are represented by (—1) in the cor- 
responding data field. The table begins with a 
header line, which attributes a descriptor to each 
column. The first column attributes a unique identi- 
fying number (ObsJD) to the observation it refers 
to. Apart from line ratios or fluxes, additional 
columns (20 to 25) are provided that describe the 
particularities of the observation (observed position, 
slit orientation, ...). The entries contained in these 
columns consist of numbers (or pointers), which refer 
to specifically coded information (see description be- 
low). This ensures that the information it contains 
can be easily accessed by any numerical program. A 
list of all the column headers and the type of infor- 
mation they contain is provided in Tabled 

As we mentioned above, the entries in some of 
the columns have a numerical format that deserves 
further clarification: "ObsJD" is the positive, con- 
secutive and unique number for row identification. 
"Position" is the value in kpc of the position within 
each galaxy. In case the observations was carried 
out perpendicular to the galaxy plane, this is the 
distance from the plane, for parallel slits this is the 
position from the galaxy center and for face-on galax- 
ies the position is determined with respect to some 
reference point that the authors had specified. The 
ratio Ha/H/3 is that reported by the authors. How- 
ever, in cases where no value was reported, a the- 
oretical ratio of 2.86 was assumed when its value 
was required to normalize the forbidden lines. The 
"EM" is the Ha emission measure, specified in units 
of pccm -6 . "ReLID" is the integer descriptor that 
identifies the bibliographical reference (see Table [lj. 
"Morphology" refers to the galaxy's morphology and 
its orientation, as follows: for spiral galaxies 11 is re- 
served for edge-on objects and 12 for face-on ones; 
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TABLE 3 

LIST OF GALAXIES CONTAINED IN DIGEDA 



Gal_ID a 


Galaxy 


R.A. 


DEC. 


]VIorphology Type 


Inclination 






(h m s) 


(d m s) 




(degrees) 


1 


M 31 


00 42 44.3 


+41 16 09 


Sb, LINER. 


77.5 


2 


M 33 


01 33 50.9 


+30 39 36 


Scd, H II 


56 


3 


M 51 


13 29 55.7 


+47 13 53 


Sbc 


64 


4 


M 81 


09 55 33.2 


+69 03 55 


Sab;LlNER Syl.8 


58 


5 


M 101 


14 03 12.6 


+54 20 57 


SBcd 


17 


20 


NGC 891 


02 22 33.4 


+42 20 57 


Sb, HII 


64 


21 


NGC 1569 


04 30 49.0 


+64 50 53 


Irr 




22 


NGC 1963 


05 32 16.8 


-36 23 55 


Scd 


85 


23 


NGC 2403 


07 36 51.4 


+65 36 09 


SBcd, H II 


62 


24 


NGC 2820 


09 21 45.6 


+64 15 29 


SBc 


90 


25 


NGC 3044 


09 53 40.9 


+01 34 47 


SBb 


84 


26 


NGC 3079 


10 01 57.8 


+55 40 47 


SBc, liner, Sy 2 


60 


27 


NGC 3628 


11 20 17.0 


+13 35 23 


Sb 


87 


28 


NGC 4013 


11 58 31.4 


+43 56 48 


Sb, HIILINER 


90 


29 


NGC 4217 


12 15 50.9 


+47 05 30 


Sb 


86 


30 


NGC 4302 


12 21 42.5 


+14 35 54 


Sc 


90 


31 


NGC 4395 


12 25 48.8 


+33 32 49 


Sm, LINER, Sy 1.8 


38 


32 


NGC 4402 


12 26 07.5 


+13 06 46 


Sb 


74 


33 


NGC 4631 


12 42 08.0 


+32 32 29 


SBb 


85 


ol 


IN V i ^tOO^ 


1 9 A 9 AC\ Q 




OC 




35 


NGC 5775 


14 53 57.6 


+03 32 40 


Sb 


86 


36 


NGC 5777 


14 51 17.8 


+58 58 41 


Sb 


83 


70 


UGC 2092 


02 36 31.6 


+07 18 34 


Scd 


86 


71 


UGC 3326 


05 39 37.1 


+77 18 45 


Scd 


90 


72 


UGC 4278 


08 13 58.9 


+45 44 32 


SBd 


90 


73 


UGC 10288 


16 14 24.8 


-00 12 27 


Sc 


87 


90 


ES0245 - G05 


01 45 03.7 


-43 35 53 


Irr 




91 


Gr 8 


12 58 40.4 


+14 13 03 


Irr 




92 


DDO 53 


08 34 07.2 


+66 10 54 


Irr 





a Unique galaxy identification number within DIGEDA. 



2 is used for irregulars. "Slit" is the slit orienta- 
tion with respect to the galactic plane: 1 for parallel 
slits, 2 for perpendicular slits and 3 for observations 
of face-on objects. "RegionJD" refers to the three 
possible emission region types: 1 and 11 are used for 
HII regions, 2 and 21 for transition zones, 3 and 31 
for DIG regions, and finally, 4 is reserved for data 
that do not fit any classification. Single digit en- 
tries denote the authors original classification (when 
it was reported) while double digit entries denote our 
own classification taking into account the value of the 
height above the galactic plane \z\ as described below 
in § 12.21 Finally, "GaUD" is an integer number that 
relates uniquely the galaxy's name to its source cata- 
log, as follows: Messier: 1-19, New General Catalog: 
20-69, Uppsala General Catalog: 70-89 and others: 
90-99. The list of all the galaxies in DIGEDA, with 
their basic characteristics, is given in Table [H 



To sum up, DIGEDA contains 1061 spectro- 
scopic observations. It is worth noting, however, that 
only 114 observations contain simultaneous measure- 
ments of [N II], [O III], [S II] and the EM (these turn 
out to be edge-on spiral galaxies). But if we consider 
instead a single pair of lines, for example [N II] and 
[O III] , the number of available data entries in this 
case increases to 462 observations. 

DIGEDA will be made available on-line, at the 
VizieR website ( |http : //vizier . u-strasbg . f rj ). 

2.2. An additional grouping criterion: \z\ 

Due to the fact that there is no standard criterion 
to distinguish DIGs from HII regions, most of the 
authors in the field report their spectroscopic studies 
without any classification. For instance, ~ 76% of 
the 1061 observations in the DIGEDA are neither 
referred to as HII region, DIG or transition zone. 
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About ~ 97% of these data correspond to edge-on 
galaxy observations, which are of particular interest 
among diffuse ionized gas studies. For this galaxy 
type, the halo gas is easier to separate from the disk 
(where most H II regions are located) . In that case, 
the DIG emission can be reasonably assumed to be 
free from contamination by star formation emission. 
Considering that in the Milky Way, the thin disk 
is defined as the region where most of the galaxy's 
star formation is taking place, we propose in this 
Paper a 'zoning' classification criterion based on the 
height (above or below) the galactic plane \z\. The 
proposed zoning classification leads to the following 
three situations: 

• \z\ < 500 pc. Owing to the fact that in the 
Milky Way, one considers that H II regions are 
confined to a disk of 350 pc height, we have 
assumed this to be a general behavior in spi- 
ral galaxies. We will therefore consider that 
line emission from regions within 500 pc from 
the galactic plane corresponds overwhelmingly 
to emission from star forming regions. Conse- 
quently, we will classify these regions as HII re- 
gions ( "Region JD" = 11). 

• 500 < \z\ < 1000 pc. The majority of spiral 
galaxies from the database do not have an incli- 
nation of exactly 90° and it is therefore prob- 
able that, due to projection effects, a line of 
sight below or above the galactic plane of be- 
tween 500 and 1000 pc will cross a large volume, 
which likely will contain both DIG and H II re- 
gion emission. Since we do not know precisely 
the contribution from each type, we will clas- 
sify these regions as being from transition zones 
( "Region JD" =21). 

• \z\ > 1 kpc. Even though some authors have re- 
ported the existence of H II regions as far away 
as 1.5 kpc from the plane, this remains an un- 
common occurrence. In this work, we will as- 
sume that the emission from any lines of sight 
above \z\ = lkpc from the galactic plane is 
free from any contamination from star forma- 
tion emission. This leads us to classify these 
regions as eDIG ("RegionJD"=31). 

Applying this 'zoning' classification to the 
database resulted in 309 HII regions, 218 transition 
zones, and 509 eDIGs. Only 2% of the observations 
could not be classified. They turn out to come from 
face-on galaxies where our criterion based on \z\ is 
not applicable. 



3. ANALYSIS AND INTERPRETATION OF THE 
DIGEDA DATA 

We have carried out a statistical study of the 
database with the aim of identifying the main dif- 
ferences between DIGs and HII regions. To assist 
us, we used the line ratio diagnostic diagrams pre- 
scribed by Baldwin, Phillips, & Terlevich (1981) and 
by Veilleux & Osterbrock (1987). Such diagrams 
were created with the intention of classifying line 
emitting regions according to their physical condi- 
tions, yet using only the most commonly measured 
lines, that is: Ha, H/3, [NIL], [Ol], [0111], and [S II] . 
Baldwin et al. (198ljf| and Veilleux & Osterbrock 
(1987) found that these diagrams provide a quantita- 
tive classification scheme of extragalactic sources, in 
which the main excitation mechanism operating on 
the gas is identified, namely: photoionization by O 
and B stars, photoionization by a power-law contin- 
uum source, shock- wave heating and photoionization 
by old and very hot stars. Apart from being used for 
the above purpose, these diagrams also turn out to 
be useful in other situations or for other kinds of 
emission objects (e.g., distinguishing planetary neb- 
ulae from HII regions, see Kniazev, Pustilnik, & 
Zucker 2008). Since the emission lines quoted above 
are among the most commonly observed lines in DIG 
regions as well, we adopted the line-ratio diagnostic 
diagrams in our statistical study. 

3.1. Results from line ratio diagnostic diagrams 

In Figures Q] and [21 we show, for the whole 
database, the line ratios [OlII]/H/3, [Ol]/Ha, 
[Sll]/Ha, and [Sll]+/Ha as a function of the ra- 
tio [Nil] /Ha. Figure [1] distinguishes the different 
galaxy types as follows: Irregular galaxies (aster- 
isks), spiral face-on galaxies (triangles) and edge-on 
spirals (squares). Figure [21 which is discussed below, 
distinguishes among the three emission region types 
defined in § O 

Our Panel (a) in Figure Q] shows that the data de- 
scribes a "seagull" shape quite similar to that found 
by Stasihska et al. (2006) from a study of 20,000 
"normal" star-forming and AGN galaxies extracted 
from the SLOAN database. While the line ratios 
from both spiral types (edge-on and face-on) extend 
over the full seagull area in Panel (a) of Figure [TJ 
those from irregular galaxies tend to occupy only 
the extreme left wing (lower [Nil] /Ha and greater 
[OlII]/H/3 values). In the other Panels (b), (c) and 
(d) of Figure [1] the same remarkable behavior is 

3 Hereafter, we will refer to the [Nll]/Ha vs. [Olirj/H/J 
diagram as 'BPT diagram'. 
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Fig. 1. Diagnostic diagrams for all observations in the database. The data points identify the galaxy types as follows: 
irregular galaxies (asterisks), face-on spiral galaxies (triangles), and edge-on spiral galaxies (squares). 



found. Notwithstanding the reasons for such behav- 
ior (see § H|) , it is statistically apparent that the data 
from irregular and spiral galaxies should not be ana- 
lyzed without differentiating the galactic type. Since 
DIGEDA contains only nine observations of irregu- 
lar galaxies, we chose to consider only spiral galaxies 
in the analysis that follows. 

Figure [2] is similar to the previous figure, except 
that it distinguishes among the three types of emis- 
sion line regions: HII regions (blue diamonds), tran- 
sition zones (green crosses) and DIGs (red pluses) . In 
Panel (a) , DIGs and H II regions conform together a 
"seagull" shape. However, it becomes evident in this 
BPT diagram that the DIGs occupy a very differ- 
ent location with respect to the HII regions. While 
H II regions are principally located in the left seagull 
wing, DIGs cover the right seagull wing, i.e. the two 
types can be distinguished from their [N II] /Ha ratio 
alone. 

Although less apparent, we can see in Panels (b) 
and (c) of Figure [2] that DIGs tend furthermore to 



have [S II] /Ha ratios greater than those from H II re- 
gions. This behavior has been previously mentioned 
on several occasions (Bland-Hawthorn et al. 1991a; 
Zurita et al. 2000; Rand 1998; Haffner, Reynolds, 
& Tufte 1999; Tullmann & Dcttmar 2000) for dif- 
ferent galaxies using different classification criteria. 
The importance of this figure in DIG studies is that 
it makes clear that even when the data have been 
classified using a very different criterion, the inten- 
sities of [N II] and [S II] are significantly greater in 
general in DIGs than in H II regions. Perhaps the 
spread in position of each type and the extent to 
which they overlap (which differ in each panel) could 
also be the manifestation of the absence of an uni- 
versally accepted classification scheme to distinguish 
the emission region type. 

Conspicuously, Panel (d) of Figure [2] does not 
show any clear trend in the behaviour of the [O I] /Ha 
line ratio. The data dispersion and the overlap in 
this panel are both larger than found with [N II] or 
[S II] , possibly because the [O I] line is the weakest 
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of all lines considered here. In this case, we can ex- 
pect the [O I] measurements to be more difficult and 
uncertain, which is consistent with the fact that we 
only have 148 [O I] measurements (all types together) 
available among the 1061 observations. 

3.2. The [N Ilj/Ha line ratio criterion 

Even though our diagnostic diagrams in Figure [2] 
strongly suggest the existence of a clear separation 
between HII regions and DIGs in the behavior of 
[Nil] /Ha and [Sll]/Ha, a significant overlap is nev- 
ertheless present. It thus becomes essential to make 
a quantitative comparison of the distributions of the 
various line ratios. This exercise would allow us to 
identify which ratios show the most significant dif- 
ferences between DIGs and H II regions in the distri- 
bution of their values. For this purpose, we can use 
box diagrams, which summarize in a very schematic 
fashion the intrinsic distribution of any of the data 
entries. In such diagrams, each box contains 50% of 
the data and each has been defined using quartiles 
Ql and Q3, which are the lowest and highest 25% 
cuts-offs in the data distribution. In such boxes, a 
vertical line denotes the position of the median of 
the distribution, while "whiskers" are used to rep- 
resent distribution boundaries that include 80% of 
the data, which have been defined using the deciles 
Dl and D9, corresponding to the lowest and highest 
10% cut-offs. 

Figure [3] shows such box diagrams for the EM 
values and for each of the most prominent line ratios 
(from spiral galaxies) contained in DIGEDA. Each 
datum's statistic distinguishes whether it considers 
only DIGs or H II regions, using the label (D) or (H), 
respectively. It is apparent that a few conspicuous 
line ratios show a distribution that is significantly 
distinct between DIGs and H II regions. These ratios 
correspond to those boxes that do not overlap (or 
overlap slightly) in log (I) position when comparing 
the (D) and (H) boxes of the same quantity. It is 
also obvious that, conversely, many line ratios show 
an important overlap in the position of the (D) and 
(H) boxes. 

If we consider the line ratios of [O I], [S II], [S II] + 
and [Nil] with respect to Ha, the following sample 
sizes are available: [O I] (157 entries), [S II] (521 en- 
tries), [SII]+ (322 entries) and [Nil] (888 entries). 
It is noteworthy that for each of these 4 line ratios, 
the distributions that the boxes in Figure [3] repre- 
sent do not overlap between the (D) and (H) entries. 
That is, these 4 line ratios show clear differences in 
their quantiles distribution between DIGs and H II 
regions. This quantitative analysis gives us further 



confidence in proposing that a few line ratios are 
by themselves sufficient to empirically distinguish a 
DIG region from an H II region, notwithstanding the 
lack of any prior knowledge about the electron den- 
sity or column size (both quantities being unknown 
in DIGs). 

Even considering that any of the above four lines 
ratios could be a useful criterion for distinguish- 
ing DIGs, the [Nil] /Ha statistic corresponds to the 
largest data sample. Observationally, it is also per- 
haps the ratio that is easiest to measure (since [N II] 
is a strong line and the [Nil] /Ha ratio does not re- 
quire any reddening correction). 

In Table [4j the numerical values of the different 
quantiles for [Nil] /Ha are sumarized. It can be con- 
cluded that emission regions with [Nil] /Ha < —0.5 
should be classified as H II regions, while those with 
[N II] /Ha > -0.3 should be classified as DIGs. Mea- 
surements of intermediate values result in an uncer- 
tain classification. In any case, we can infer that 
emission regions with a ratio [Nil] /Ha < —0.4 have 
a high probability of being an H II region while re- 
gions with [Nil] /Ha > —0.4 have a higher probabil- 
ity of being a DIG. 

Additionally, the boxes in Figure [3] representing 
the ratios [OlII]/H/3, [NII]/[SII], [NII]/[SII] + , and 
[NII]/[OlII] do not show a significant difference in 
distribution between DIGs and HII regions. These 
ratios are therefore useless for distinguishing one 
type of region from another. 

Some authors have proposed to use the EM as 
a classification criterion (e.g., Walterbos & Braun 
1994). We note, however, from Figure[3]that the box 
diagrams of EM do not show any significant differ- 
ences between the DIG and H II region distributions. 
Other authors (e.g., Bland-Hawthorn et al. 1991a, b; 
Rand 1998) have suggested that the [Nll]/Ha ra- 
tio shows a strong anti-correlation with the EM, i.e. 
the [Nil] /Ha would increase as the EM decreases. 
To which extent therefore can the EM be a distin- 
guishing universal criterion for DIGs? Even though 
the EM distribution box in Figure[3]suggests that we 
cannot (because of a large overlap of the (D) and (H) 
boxes), we may nevertheless consider the possibility 
that within a single galaxy, this criterion might be 
applicable. After all, it is the only quantity in the 
box diagrams of Figure [3] that does not correspond 
to a ratio and is therefore not normalized. The wide 
overlap might for instance be the result of a wide 
spread in EM values from galaxy to galaxy. 

In order to explore further this question, we plot- 
ted EM values against several emission line ratios. 
As expected, we found a very feeble general rela- 
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tionship and only in the case of [N Il]/Ha. However, 
when we considered individual galaxies separately, a 
different picture emerged. Our results are summa- 
rized in Figure |H where Panel (a) has [Nil] /Ha in 
the abscissa and Panel (b) [OlII]/H/3, while the EM 
is displayed in the ordinate. A different color cod- 
ing is used for each galaxy: M33 (magenta), M51 
(cyan), NGC 1963 (yellow), NGC 3044 (red), NGC 
4402 (green), NGC 891 (black), NGC 4634 (blue, 
Panel (a) only), NGC 4631 (blue, Panel (b) only). 
The best log-log fits of the points of each galaxy are 
overplotted. In Panel (a), we see that there exists 
indeed a relationship between EM and [N II] /Ha for 
each galaxy individually. In the case of [O III] /H/3 
(Panel b) , a similar but weaker relationship appears 
to be present as well. 

Even though a general EM-based criterion cannot 
be defined (Figure |3J , it is apparent from Figure 0] 
that when the [Nil] /Ha ratio (Table 0} is used to 
determine the emission region type, a critical value 
EM C can then be found that appears to be different 
for each galaxy individually. We obtain the interest- 
ing result that emission regions with EM > EM C can 
be considered to be HII regions while emission re- 
gions in the same galaxy with EM < EM C turn out to 
be DIGs. This particular EM C criterion, defined for 
each galaxy, is based on the more general [N II] /Ha 
criterion. EM C allows us now to distinguish DIGs 
from H II regions with certainty, confirming that the 
basic physical conditions (as represented by line ra- 
tios) in these two region types are indeed different. 

4. DISCUSSION 

Following our statistical analysis of the database, 
there are two results that deserve further discussion: 
on the one hand, the different physical conditions 
taking place in DIGs and H II regions (in the con- 
text of spiral galaxies) and, on the other hand, the 
extreme position that the DIGs from irregular galax- 
ies occupy in the BPT-diagram with respect to the 
DIGs from spirals. 

4.1. DIG and H II regions: what is the underlying 
difference? 

BPT proposed a way to classify emission line 
regions, which is based on their location in the 
[O III]/H/3 vs. [N II] /Ha diagnostic diagram. These 
authors were successful in separating H II regions 
(low [N II] /Ha) from liners and Seyfert galaxies 
(high [Nil] /Ha). Using data from DIGEDA, we 
found that most of the DIGs fall in the regions oc- 
cupied by liners and Seyfert s (see Figure How- 
ever, we cannot foresee how in non-active galaxies 



the same kind of photoionization source that pow- 
ers AGN could replace O-B stars of H II regions and 
become responsible for the excitation of the DIG. 

More recently, using the SDSS survey Stasihska 
et al. (2008) proposed an original interpetation for 
objects occupying the right part of the BPT dia- 
gram (which they dubbed the seagull's right wing). 
They pointed out that the simultaneous increase of 
[OlII]/H/3 and [Nil] /Ha can also be the result of 
photoionization by hot old stars, for example the 
central stars of planetary nebulae (PNe). According 
to this interpretation, the ionizing source of the DIG 
would be similar to that of H II regions (i.e. a stel- 
lar continuum), except for the effective temperature 
of the ionizing stars and, consequently, the higher 
electron temperature of the gas. 

Following the conclusions reached by Stasihska et 
al. (2008), we propose that the DIG could be ion- 
ized by old stars, at least in spiral galaxies. In fact, 
the suggestion that a hot stellar population might 
contribute to the ionizing source was first proposed 
by Lyon (1975) and later explored by Sokolowski & 
Bland-Hawthorn (1991) who developed photoioniza- 
tion models that combine OB stars with old hot 
stars. Sokolowski & Bland-Hawthorn (1991) con- 
cluded that such composite models can reproduce 
the DIGs observations in NGC 891. However, as 
early as in 1991, the available observational data 
were scant and no [OlII]/H/3 ratios from DIGs had 
yet been published. 

One can find other explanations for the position 
of DIGs in the BPT diagram besides that of a hotter 
photoionization source. For instance, the photoion- 
ization models developed by Sokolowski (1993) sug- 
gest that the ionizing flux from O and B stars could 
simulate the effect of a hotter source when interven- 
ing dust is present (between the stars and the DIG) 
and the selective extinction by such dust is taken 
into account. A similar UV hardening could also 
result from the radiation that leaks out from disk 
HII regions (Domgorgen & Mathis 1994; Zurita et 
al. 2000). Both effects of selective absortion and of 
photoionization by old and hot stars are the subject 
of a paper in preparation. 

4.2. DIG in irregular galaxies 

The DIGs observed in irregulars occupy the same 
part of the BPT-diagram as do the H II regions of 
spiral galaxies (low [Nil] /Ha in Figure [T]). Fur- 
thermore, we have no evidence of any DIG from ir- 
regulars occupying the high [N II] /Ha part of the 
diagram. About the first aspect mentioned above 
(low [Nil] /Ha), given the fact that the position of 
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Fig. 2. Diagnostic diagrams for all observations in the database. The data points identify the emission region types as 
follows: H II regions (blue diamonds), transition zones (green crosses) and DIGs (red pluses). 



an emission line object in the BPT diagram is re- 
lated to specific physical conditions of the gas, our 
results suggest that the so-called DIG from irregular 
galaxies could in fact just be the manifestation of 
low surface brightness H II regions. 

On the other hand, the lack of DIG from irreg- 
ulars at high [Nil] /Ha could have several alterna- 
tive explanations. First, it could be that this emis- 
sion is too weak to have been detected yet. Another 
possible explanation follows from the idea developed 
in the previous section, namely, that the DIG ion- 
ization results from "old stars". It is conceivable 
that in irregular galaxies the old stellar population 
is proportionally much smaller than in spiral galax- 
ies, the effect of which would be that this population 
could not then contribute significantly to the ioniza- 
tion budget. In any event, since we only have a few 
data points in the case of irregular galaxies (7 H II 
regions and only 2 DIGs with measurements of [N II] 
and [O III] ) , we cannot yet reach definitive conclu- 
sions about irregulars. 



TABLE 4 
QUANTIL VALUES FOR [Nil] /Ha 





Dl 


Ql 


Median 


Q3 


D9 


H II reg 
DIG 


-0.89 
-0.51 


-0.68 
-0.38 


-0.52 
-0.22 


-0.41 
-0.04 


-0.33 
0.09 



5. CONCLUSIONS 

In this work, we present the first emission line 
database from DIGs (DIGEDA) made up from 
spectroscopic data available in the literature. The 
database is a compilation of 17 bibliographical ref- 
erences. It contains 1061 observed regions (309 HII 
regions, 218 transition zones and 509 DIGs) out of 
29 galaxies (spirals and irregulars) and is freely avail- 
able from CDS. 

DIGEDA has allowed us to carry out for the first 
time a statistical analysis that aimed at character- 
izing the general behavior of the strongest emission 
lines, and at finding global trends among the vari- 
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Fig. 3. Box diagrams of the data distributions of the EM 
and of the most prominent line ratios from DIGs and H II 
regions (D and H labels, respectively) . All the values are 
logarithmic, with that of the EM reduced by 3 dex. 



ous line ratios, including the EM. One of the results 
is that the DIGs of irregular galaxies show extreme 
values in their line ratios with respect to that of spi- 
rals (see §[3] and §2]). These, in fact, lie close to the 
values observed in H II regions. Since the available 
data are very limited, it is indispensable to obtain 
more observations to confirm this result. 

The analysis carried out with DIGEDA leads 
us to define a universal criterion for spiral galax- 
ies, which allows distinguishing DIGs from H II re- 
gions. The distribution of [Nil] /Ha ratios shows 
that [Nil] /Ha = —0.4 defines a critical value: any 
emission region with a ratio greater than —0.3 should 
be classified as DIG while a ratio below —0.5 is likely 
to come from an H II region. 

Taking advantage of the high number of data in 
DIGEDA, we could confirm (or refute) the reality of 
line ratio behaviors that have been proposed in the 
literature concerning individual galaxies. We could 
check to what extent these applied to all DIGs taken 
together. This is the case for the [Nil] /Ha ratio, 
which shows a well defined anti-correlation with EM, 
but of varying slopes between galaxies. These varia- 
tions in the log-log slopes rule out any definition of 
a universal EM value to distinguish DIGs from H II 
regions. Nevertheless, a critical value, EM C , can be 
defined for each galaxy individualy. 
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Fig. 4. The EM against (a) [Nil] /Ha and (b) [Olll}//3. 
A different color coding is used for each galaxy: M33 
(magenta), M51 (cyan), NGC 1963 (yellow), NGC 3044 
(red), NGC 4402 (green), NGC 891 (black) and NGC 
4634 (blue, in panel a), NGC 4631 (blue, in panel b). 
The best log-log linear fits of the points of each galaxy 
are overplotted using the same color coding. 
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